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ABSTRACT 

We report the results of a multi-band observing campaign on the famous blazar 3C 279 conducted during a 
phase of increased activity from 2013 December to 2014 April, including first observations of it with NuSTAR. 

The 7 -ray emission of the source measured by Fermi-LACT showed multiple distinct flares reaching the highest 
flux level measured in this object since the beginning of the Fermi mission, with F{E > lOOMeV) of 10“^ 
photons cm“^ s~\ and with a flux doubling time scale as short as 2 hours. The 7 -ray spectrum during one of the 
flares was very hard, with an index of = 1.7 ±0.1, which is rarely seen in flat spectrum radio quasars. The 
lack of concurrent optical variability implies a very high Compton dominance parameter L^/L^yn > 300. Two 
1-day NuSTAR observations with accompanying Swift pointings were separated by 2 weeks, probing different 
levels of source activity. While the 0.5-70keV X-ray spectrum obtained during the first pointing, and fitted 
jointly with Swift-XKT is well-described by a simple power law, the second joint observation showed an unusual 
spectral structure: the spectrum softens by AFx — 0.4 at ^ 4keV. Modeling the broad-band SED during this 
flare with the standard synchrotron plus inverse Compton model requires: (1) the location of the 7 -ray emitting 
region is comparable with the broad line region radius, ( 2 ) a very hard electron energy distribution index p ~ 1 , 

(3) total jet power significantly exceeding the accretion disk luminosity Lj/Ld > 10, and (4) extremely low jet 
magnetization with Lb/Fj < 10“^. We also find that single-zone models that match the observed 7 -ray and 
optical spectra cannot satisfactorily explain the production of X-ray emission. 

Subject headings: galaxies: active — galaxies: jets — gamma rays: galaxies — quasars: individual (3C 279) 

— radiation mechanisms: non-thermal—X-rays: galaxies 
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1 . INTRODUCTION 

Blazars are active galaxies where the strong, non-thermal 
electromagnetic emission, generally detected in all observable 
bands from the radio to 7 -ray spectral regimes, is dominated 
by the relativistic jet pointing close to our line of sight. De¬ 
tailed studies of blazar spectra, and in particular the spectral 
variability, are indispensable tools to determine the physical 
processes responsible for the emission from the jet, leading 
to understanding the distribution of radiating particles, and 
eventually, the processes responsible for their acceleration. 

3C 279 is among the best studied blazars; it is detected in 
all accessible spectral bands, revealing highly variable emis¬ 
sion. It consistently shows strong 7 -ray emission, already 
clearly detected with the EGRET instrument on the Comp¬ 
ton Gamma-Ray Observatory (CGRO; Hartman et al. 1992). 
The object, at z = 0.536 (Lynds et al. 1965), is associated with 
a luminous flat-spectrum radio quasar (ESRQ) with promi¬ 
nent broad emission lines. Optical and UV observations in 
the low-flux state (Plan et al. 1999) allow the luminosity of 
the accretion disk to be estimated at Ld ~ 6 X 10"^^ erg s 
The estimates of the mass of the central black hole are in 
the range of (3-8) x IO^Mq, derived from the luminosity 
of broad optical emission lines (Woo & Urry 2002), the width 
of the Hp line (Gu et al. 2001), and the luminosity of the host 
galaxy (Nilsson et al. 2009). The object possesses a compact, 
milliarcsecond-scale radio core and a jet with time-variable 
structure. Multi-epoch radio observations conducted between 
1998 and 2001 by Jorstad et al. (2004, 2005) provided an esti¬ 
mate of the bulk Lorentz factor of the radio-emitting material, 
rj = 15.5±2.5 and the direction of motion to the line of sight, 
0obs = 2.1 ± 1.1 degrees, which corresponds to a Dopplerfac- 
torDof24.1±6.5. 

As is the case for blazars, the most compelling mechanism 
for the production of the radio through optical bands is syn¬ 
chrotron emission, while the 7 -rays arise via inverse Comp¬ 
ton emission by the same relativistic electrons producing the 
synchrotron emission (Sikora et al. 2009). Alternative models 
involving hadronic interactions require significantly higher jet 
powers due to their lower radiative efficiency (Bottcher et al. 
2009). Since in the co-moving frame of the relativistic jet 
the photon energy density in luminous blazars is dominated 
by external radiation sources, production of 7 -rays is most 
efficient by scattering of the external photons (Dermer et al. 
1992; Sikora et al. 1994). 3C 279 is regularly monitored by 
the Fermi satellite together with many different facilities cov¬ 
ering a range of spectral bands, from radio and optical to 
X-rays. The correlations of the highly variable time series 
between the optical polarization level/angle and 7 -rays pro¬ 
vide strong evidence for the synchrotron ■¥ Compton model, 
and suggest among the solutions that the jet structure is not 
axisymmetric (Abdo et al. 2010a), or the presence of a heli¬ 
cal magnetic held component (Zhang et al. 2015). The rapid 
variability, together with the rate of change of the polariza¬ 
tion angle, suggest a compact (light days) emission region 
that is located at an appreciable (> a parsec) distance along 
the jet from the black hole. Eurthermore, the close but not ex¬ 
act correlation of the optical and 7 -ray flares, with the optical 
lagging the 7 -rays by ~ 10 days (Hayashida et al. 2012), has 
supported this basic scenario (Janiak et al. 2012). 

Perhaps the largest mystery in 3C 279 — and other lumi- 

Pian et al. (1999) report a lower value Lj = 2 X 10'*^ ergs^*, appai'- 
ently without a bolometric con'ection, and we used that value previously in 
Hayashida et al. (2012). Here, we apply a correction by factor 3. 


nous blazars as well — is the nature of its X-ray emission 
(Sikora et al. 2013). Early comparison of the RATE X-ray and 
EGRET 7 -ray time series revealed a close association of the 
7 -ray and X-ray flares (Wehrle et al. 1998), suggesting that 
the X-ray flux might be the low-energy end of the same in¬ 
verse Compton emission component detected at higher ener¬ 
gies by EGRET. This is supported indirectly by a good overall 
correlation between long term RXTE and optical data (which, 
according to the above, should be a reasonable proxy for 7 -ray 
flux), although individual flares show time lags up to ~ ±20 
days (Chatterjee et al. 2008). However, better sampling pro¬ 
vided by the multi-band time series covering many years (and 
owing mainly to the all-sky monitoring capability of the Fermi 
Large Area Telescope, EAT; Atwood et al. 2009) revealed that 
the 7 -ray and X-ray fluxes are often not well correlated both 
for this object (Hayashida et al. 2012) and for other blazars 
(see e.g., Bonning et al. 2009). The nature of blazar X-ray 
emission is still somewhat unclear. 

3C 279 is also a prominent hard X-ray and soft 
7 -ray source, detected by CGRG/OSSE (Hartman et al. 
1996), INTEGRAL (Beckmann et al. 2006), and 5w//f-BAT 
(Tueller et al. 2010). However, these observations did not pro¬ 
vide a precise measurement of the hard X-ray spectrum of 
this source, that would allow discrimination between alterna¬ 
tive spectral components. 3C 279 was selected as one of a 
few blazar targets to be observed in the early phase of the Nu¬ 
clear Spectroscopic Telescope Array (NuSTAR', Harrison et al. 
2013) focusing hard X-ray (3-79 keV) mission. 

After a brief hiatus, 3C 279 became very active late in 2013, 
producing a series of 7 -ray flares and reaching the highest 7 - 
ray flux level recorded by the Eermi-LAT (Buson 2013) for 
this source. The flaring activities of 3C 279 triggered many 
observations, including the first two pointings by the NuS¬ 
TAR satellite, enabling sensitive spectral measurements up to 
70keV. Here, we present the results of the analysis of the 
Fermi-LPsA, NuSTAR, and Swift data together with optical ob¬ 
servations by the SMARTS and Kanata telescopes as well as 
the sub-mm data from the Submillimeter array (SMA) col¬ 
lected for five months during the high-activity period, from 
2013 November to 2014 April. A part of this period — 2014 
March-April — was studied independently by Paliya et al. 
(2015). In Section 2 we describe in detail the data analysis 
procedures and basic observational findings. In Section 3 we 
compare the observational results between multiple bands. In 
Section 4 we discuss the interpretation and theoretical impli¬ 
cations of our results, and we conclude in Section 5. 

2 . OBSERVATIONS AND DATA REDUCTION 
2.1. Fermi-LAT: Gamma-ray observations 

The EAT is a pair-production telescope onboard the Fermi 
satellite with large effective area (~ 6500 cm^ on axis 
for 1 GeV photons) and a large held of view (2.4 sr), 
sensitive from 20MeV to 300 GeV (Atwood etal. 2009). 
Here, we analyzed EAT data for the sky region includ¬ 
ing 3C 279 following the standard procedure^®, using the 
EAT analysis software ScienceTools v9r34vl with the 
P7REP_SOURCE_V15 instrument response functions. The 
azimuthal dependence of the effective area was taken into ac¬ 
count for analysis with short time scales (<1 day). Events 
in the energy range 0.1-300GeV were extracted within a 15° 
acceptance cone of the Region of Interest (ROI) on the loca- 

http; / / f ermi . gsf c. nasa . gov/ssc/data/analysis/ 
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Fig. 1.— Light curves of 3C 279 in the 7 -ray band (integral photon flux) as observed by Fermi-LPC^. Top panel shows the long-term light curve above 100 MeV 
in 3-day bins. The other panels show light curves for the 2013-2014 active period: from the top to bottom, (1) above 100 MeV in 6 -hour bins, (2) from 100 MeV 
to 1 GeV in 1-day bins, (3) above 1 GeV in 1-day bins, (4) arrival time distribution of photons with energies above lOGeV, (5) photon index of 3C 279 above 
100 MeV in 1-day bins. A gap in the data around ~MJD 56680-56690 is due to a ToO observation of the Crab Nebula, during which time no exposure was 
available in the direction of 3C 279. The vertical bars in data points represent 1 cr statistical en'ors and the down aiTows indicate 95% confidence level upper 
limits. 


tion of 3C 279 (RA = 195.047°, DEC=-5.789°, J2000). It is 
known that the Sun comes very close to and occults 3C 279 
on October 8 each year. The data when the source is within 
5° of the Sun were excluded. Gamma-ray fluxes and spec¬ 
tra were determined by an unbinned maximum likelihood fit 
with gt 1 ike. We examined the significance of the g-ray sig¬ 
nal from the sources by means of the test statistic (TS) based 
on the likelihood ratio test^'. The background model included 

TS = 25 with 2 degrees of freedom corresponds to an estimated ~ 4.6 tr 
pre-trials statistical significance assuming that the null-hypothesis TS distri¬ 
bution follows a distribution (see Mattox et al. 1996). 


all known g-ray sources within the ROI from the 2-year LAT 
catalog (2FGL: Nolan et al. 2012). Additionally, the model 
included the isotropic and Galactic diffuse emission compo¬ 
nents^^. Flux normalizations for the diffuse and background 
sources were left free in the fitting procedure. 

Several g-ray light curves as measured by Fermf-LAT can 
be seen in Figure 1 . The top stand-alone panel shows the g- 
ray flux above 100 MeV for about 6 years since the begin¬ 
ning of scientific operations of the Fermi-LAT (2008 August 
5) up to 2014 August 31 (MJD 54683-56900) binned into 

iso_source_v05 .txt and gll_iem_v05_revl. fit 
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Fig. 2.— Gamma-ray light curves (integral photon flux) of 3C279 around the three large flares with fine time bins. Top panels: > 100MeV; lower panel: > 
1 GeV. For Flares 1 and 2, the bins are equal to two Fermi orbital periods (192 min). For Flare 3, during a ToO observation, the bins are equal to one Fermi orbital 
period (96 min). The vertical bars in data points represent 1 cr statistical errors and the down arrows indicate 95% confidence level upper limits. 


TABLE 1 

Fitting results of the light curve profile in the 7-ray band measured by Ferm/-LAT. 


Flare 

number 

Tise 

(hi-s) 

(firs) 

b 

(10“^ photons cm“^ s“^) 

Fo 

(10 ^ photons cm ^ s 

^0 

(MJD) 

Flare 1 

1.4±0.8 

7.4 ±3.2 

150 ±36 

19±12 

56646.35 ±0.04 

Flare 2 

6.4 ±2.4 

0.68 ±0.59 

100 ±26 

19±5 

56718.32 ±0.07 

Flare 3 (ToO) 

2 .6±0.6 

5.0±0.8 

216±19 

10.5 ±6.6 

56750.30 ±0.04 


3-day intervals. After ~MJD 56600, the source entered the 
most active state since the launch of Fermi satellite. This re¬ 
sulted in Target-of-Opportunity (ToO) pointing observations 
for 3C 279, which were performed between 2014 March 
31 21:59:47 UTC (MJD 56747.91652) and 2014 April 04 
12:42:01 UTC (MJD 56751.52918), and those observations 
are included in our analysis. The time series of the 7 -ray flux 
and photon index of 3C 279 measured with Uermi-LAT during 
the most active states from MJD 56615 (2013 November 19) 
to MJD 56775 (2014 April 28), are illustrated in other panels 
in Figure 1 . 

Three distinct flaring intervals are evident in the 7 -ray light 
curve: Flare 1 (-MJD 56650), Flare 2 (-MJD 56720) and 
Flare 3 (~MJD 56750). The maximum 1-day averaged flux 
above 100MeV reached (62.2 ±2.4) x 10“^photons cm“^ s“* 
(TS = 3892) on MJD 56749 (2014 April 03) which is 
about three times higher than the maximum 1 -day averaged 
flux recorded during the first two years (on MJD 54800: 
Hayashida et al. 2012). On the other hand, the maximum 1- 
day averaged flux above 1 GeV was observed on MJD 56645 
(2013 December 20) at (9.8 ± 1.2) x 10“^ photons cm“^ s“', 
much higher than the > 1 GeV flux on MJD 56749, which was 
(3.9±0.4) X 10“^photons cm“^ s“'. The photon index also 
shows a hardening trend toward MJD 56645, when it reached 
a very hard index of 1.82 ± 0.06, which is rarely observed in 
FSRQs. 

Figure 2 shows detailed light curves around the flares with 
short time bins. During Flares 1 and 2, the fluxes were de¬ 
rived with an interval of 192 min, corresponding to two or¬ 
bital periods of Fermi-LAT. During Flare 3, because the ToO 
pointing to 3C 279 increased the exposure, time bins as short 
as one orbital period (96 min) were used. The peak flux above 
100 MeV in those time intervals (192 min and 96 min) reached 

throughout this paper, each error represents ale statistical uncertainty 


^ 120 X 10“^ photons cm“^ s“* . 

The very rapid variability apparent in the data can be fitted 
by the following function to characterize the time profiles of 
the source flux variations: 

b 

F{t) = Fa + ——-7 -^— (1) 

g-h-tol/r-risc +g(f-'o)/Tfall 

This formula has also been used in variability studies of 
other LAT-detected bright blazars to characterize the temporal 
structure of 7 -ray light curves (Abdo et al. 2010c). The dou¬ 
ble exponential form has been applied previously to the light 
curves of blazars (Valtaoja et al. 1999) as well as gamma-ray 
bursts (e.g., Norris et al. 2000). In this function, each Trise and 
Tfaii represents the "characteristic" time scale for the rising and 
falling parts of the light curve, respectively, and fo describes 
approximatively the time of the peak (it corresponds to the 
actual maximum only for symmetric flares). In general, the 
time of the maximum of a flare {tp) can be described using 
parameters in Equation 1 as: 


tp 


AiseTfall 

- to-\ - 

Aise “t" Aall 


In 



( 2 ) 


The parameters of the fitting results are summarized in Ta¬ 
ble 1. The time profiles show asymmetric structures in all 
flares; generally the rise times correspond to 1-2 hours, which 
are several times shorter than the fall times of 5-8 hours 
in Flares 1 and 3. On the other hand, the fall time ap¬ 
pears to be less than 1 hour in Flare 2 (although the fit¬ 
ting error of the parameter is quite large). One can see in 
the light curve of Flare 2 in Figure 2 that the flux reached 

90 X 1 0“’ photons cm“^ s“' at the peak but suddenly dropped 
by a factor of ^ 3 in the next bin, 2 orbits (196 min) later. 

Gamma-ray spectra were extracted from the following four 
periods: 


• (A) Overlapping with the first NiiSTAR observation 
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TABLE 2 

Results of spectral fitting in the 7-ray band measured by E^rm/-LAT. 


Period 

(MJD - 56000) 

fitting modeP 

Gainma-ray spectnam {Fermi-LAl) 

r^/a/r^i ^/r^2 £^brk(GeV) 

TS 

-2AA 

Flux (> 0.1 GeV) 
( 10 ^^ photons crrT^ s“'; 

# ot photons 
> 10 GeV 

Penod A (3 days) 

PL 

2.36±0.13 



174 


5.9 ±0.9 

1 

Dec. 16, Oh - 19, Oh 
(642.0-645.0) 

LogP 

2.32±0.17 

0.03 ±0.07 


174 

< 0.1 

5.7 ±0.9 

(26.1 GeV) 

Period B (0.2 days) 

PL 

1.71±0.10 





117.6± 19.7 

1 

Dec. 20,9h36 - 14h24 

LogP 

1.12±0.31 

0.19 ±0.09 


413 

6.0 

94.5 ±18.1 

(10.4 GeV) 

(646.4—646.6) 

BPL 

1.41 ±0.17 

3.01 ±0.91 

3.6±1.6 

415 

7.6 

100.6±18.4 


Period C (3 days) 

PL 

2.29±0.13 



219 


17.1±2.8 

1 

Dec. 31, Oh - Jan. 03, Oh 

LogP 

2.29±0.16 

0.00 ±0.06 


219 

< 0.1 

17.1±2.9 

(14.7 GeV) 

(657.0-660.0) 

BPL 

2.22 ±0.42 

2.32±0.20 

0.34 ±0.27 

219 

< 0.1 

16.9±3.1 


Period D (0.267 days) 

PL 

2.16±0.06 



1839 


117.9±7.1 

1 

Apr. 03,5h03 - 1 lh27 

LogP 

2.02 ±0.08 

0.10±0.05 


1840 

5.3 

114.9±7.1 

(13.5 GeV) 

(750.210-750.477) 

BPL 

2.02 ±0.09 

2.89 ±0.45 

1 . 6 ± 0.6 

1843 

8.0 

115.1±7.7 



^ PL: power law model, LogP: log parabola model, BPL: broken power law model. See definitions in the text. 

AL represents the difference of the logarithm of the total likelihood of the fit with respect to the case with a PL for the source. 



Fig. 3.— Gamma-ray spectral energy distribution of 3C279 as measured 
by Fermi-LP^ during the four periods identified in the text (see Section 2.1) 
as well as in Table 2. The plot includes the spectra of 3C 279 from the 2008- 
2010 campaign (Hayashida et al. 2012), including a large flare and a two-year 
average. In data points, the horizontal bars describe the energy ranges of bins 
and the vertical bars represent 1 a statistical errors. The down ari'ows indicate 
95% confidence level upper limits. 

(see Section 2.2.1). Although the NuSTAR observation 
lasted for about one day, in order to increase the 7 -ray 
photon statistics, the LAT spectrum was extracted from 
3 days where the source showed comparable flux level 
(as inferred from the light curve with 1-day bins). In 
this period, the source was found to be in a relatively 
low state. 

• (B) For three orbits (~ 4.5 hours) at the peak of Flare 1, 
when the source showed a very hard 7 -ray photon index 
« 2 ). 

• (C) Overlapping with the second NuSTAR observation 
(see Section 2.2.1). As in the case of Period A, the 
length of this period is 3 days, while the NuSTAR obser¬ 
vation lasted about 1 day. The source flux was higher 
than in Period A. 

• (D) At the peak of Flare 3 for 4 orbits (^6 hours). 


In a similar manner to previous spectral studies of 
the source with the Ferm/-LAT (Hayashida et al. 2012; 
Aleksic et al. 2014a), each 7 -ray spectrum was modeled us¬ 
ing a simple power-law (PL; dN/dE oc a broken 

power-law (BPL; chV /dE oc for E < Lbrk and chV /dE oc 

otherwise), and a log-parabola model (LogP; chV /dE (x 
with Eq = 300 MeV). The spectral fit¬ 
ting results are summarized in Table 2 and a spectral en¬ 
ergy distribution for each period is plotted in Figure 3. In 
contrast to the general feature of FSRQs that the photon 
index is almost constant regardless of the source flux (see 
e.g., Hayashida et al. 2012), the spectral shape significantly 
changed between the periods. Remarkably, the photon in¬ 
dex of the simple PL model for the Period B resulted in 
an unusually hard index for FSRQs, of 1.71 ±0.10. Such 
a hard photon index has not been previously reported in 
past LAT observations of 3C 279 that included several flar¬ 
ing episodes (Hayashida et al. 2012; Aleksic et al. 2014a). 
Among the sources in the Second LAT AGN Catalog (2LAC: 
Abdo et al. 2010b), the mean photon index value of FSRQs is 
2.4, and only one FSRQ (2FGL J0808.2-0750) in the Clean 
and flux-limited Sample has the photon index of F.^ <2 (see 
Figure 18 in Abdoet al. 2010b). Occasionally, hard photon 
indices have been observed in bright FSRQs during rapid flar¬ 
ing events (Pacciani et al. 2014). The photon index of Period 
B is even harder than the index of 4CH-2L35 (1.95 ±0.21; 
Aleksic etal. 2011b) and of PKS 1510-089 (2.29 ± 0.02; 
Aleksic et al. 2014b) at the time when the > 100 GeV emis¬ 
sion was detected. 

No significant deviations from a PL model were detected in 
the spectra of Period A and C, while evidence of spectral cur¬ 
vature was observed in the spectra of the flare peaks. Period 
B and D. As derived fitting the BPL model for Period B, the 
photon index of the lower energy part (below Lbrk = 3.6 ± 1 .6 
GeV) is 1.41 ±0.17. This is comparable to the photon index 
of the rising part of the inverse-Compton emission for the case 
of a parent electron index of 2 , as is typical of the 7 -ray spec¬ 
tra of high-frequency peaked BL Lac objects. One can easily 
recognize such a rising spectral feature in Figure 3. On the 
other hand. Period D also shows a very high flux, exceeding 
10 “^ photons cm“^ s“'(> lOOMeV), comparable to the flux of 
Period B. However, spectral shape is characterized by a soft 
index (F.^ > 2). The photon index of the lower energy part as 
derived from fitting with the BPL model is not significantly 
harder than 2 , nor is a rising spectral feature apparent in the 
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SED plot of Figure 3. 

2.2. X-ray observations 
2.2.1. NuSTAR: hard X-rays 

NuSTAR is a Small Explorer satellite sensitive to hard 
X-rays, covering the bandpass of 3 - 79keV. It features 
two multilayer-coated optics, focusing the reflected X-rays 
onto CdZnTe pixel detectors which provide spectral resolu¬ 
tion (FWHM) of 0.4 keV at lOkeV, increasing to 0.9 keV at 
68keV. The held of view of each telescope is ^ 13', and the 
half-power diameter of the point spread function is ^ 1'. The 
low background resulting from focusing of X-rays provides 
an unprecedented sensitivity for measuring fluxes and spec¬ 
tra of celestial sources. For more details, see Harrison et al. 
(2013). 

NuSTAR observed 3C 279 twice. The first observation was 
performed between 2013 December 16, 05:51:07 and 2013 
December 17, 04:06:07 (UTC) and the second one between 
2013 December 31, 23:46:07 and 2014 January 01, 22:11:07 
(UTC). The raw data products were processed with the NuS¬ 
TAR Data Analysis Software (NuSTARDAS) package v.1.3.1, 
using the nupipeline software module which produces 
calibrated and cleaned event files. We used the calibration 
files available in the NuSTAR CALDB calibration data base 
V.20140414. Source and background data were extracted from 
a region of U5 radius, centered respectively on the centroid 
of the X-ray source, and a region 5' N of the source loca¬ 
tion on the same chip. Spectra were binned in order to over¬ 
sample the instrumental resolution by at least a factor of 2.5 
and to have a signal-to-noise ratio greater than 4 in each spec¬ 
tral channel. Net “on source” exposure times corresponded 
to ~ 39.6ks and ~ 42.7ks for the first (Dec. 16) and second 
(Dec. 31) observations, respectively. We considered the spec¬ 
tral channels corresponding nominally to the 3.0-70 keV en¬ 
ergy range. The net (background subtracted) count rates for 
the first observation were 0.303 ±0.003 and 0.294 ± 0.003 cnt 
s“^ respectively for module A and module B, while for the sec¬ 
ond observation they were 0.636 ±0.004 and 0.590 ±0.004 
cnt s“*. We plotted the raw (not background-subtracted) 
counts binned on an orbital time scale in Figure 4. It is ap¬ 
parent that the source was variable from one observation to 
the other, but also that the source varied within the second 
observation via secular decrease of flux for the first 3 hours, 
followed by an increase by nearly a factor of 2. 

The spectral fitting was performed using XSPEC vl2.8.1 
with the standard instrumental response matrices and effec¬ 
tive area files derived using the NuSTARDAS software mod¬ 
ule nuproducts. For each observation, we fitted the two 
modules simultaneously including a small normalization fac¬ 
tor for module B with respect to the module A in the model pa¬ 
rameters. We adopted simple power-law and a broken-power- 
law models modified by the effects of the Galactic absorption, 
corresponding to a column of 2.2 x 10^° cm“^ (Kalberla et al. 
2005). The results of the two spectral fits were compared 
against each other by using an F-test to examine improve¬ 
ments by the broken-power-law model. The simple power- 
law model gave acceptable results for both observed spec¬ 
tra, with x^/d o.f. of 666.8/660 (41.9% for the correspond¬ 
ing probability) and 831.2/886 (90.6%), respectively. Al¬ 
though the model fluxes in the 2-10 keV band between the 
two observations showed a difference of about a factor of 2, 
the resulting photon indices were similar: 1.739 ± 0.013, and 
1.754 ± 0.008. While there was no improvement in the fit ob¬ 


tained using the broken-power-law model in the first observa¬ 
tion, it gave slightly better fits for the spectrum of the second 
observation, with a x^/d.o.f. of 820.8/884 (93.6%), yielding 
a probability of 0.4% (~ 2.9 a) that the improvement in the 
fit was due to chance (as assessed with an F-test). This may 
indicate a deviation from a single power law in the spectrum 
of the second observation. The fitting results for the NuSTAR 
spectra are summarized in Table 3. 

2.2.2. Swift-XRT: X-ray 

The publicly available Swift X-Rny Telescope (XRT) data in 
the HEASARC database^'^ reveal that Swift observed 3C 279 
51 times between 2013 November and 2014 April. We an¬ 
alyzed all those observation IDs (ObsIDs). The exposure 
times ranged from 265 s (ObsID:35019120) to 9470 s (Ob- 
sID:35019100). The XRT was used in photon counting mode, 
and no evidence of pile-up was found. The XRT data were 
first calibrated and cleaned with standard Altering criteria 
with the xrtpipeline software module distributed with 
the XRT Data Analysis Software (XRTDAS, version 2.9.2). 

The calibration files available in the version 20140709 of the 
Swift-XRT CAFDB were used in the data reduction. The 
source events were extracted from a circular region, 20 pixels 
(1 pixel ~ 2."36) in radius, centered on the source position. 

The background was determined using data extracted from a 
circular region, 40 pixels in radius, centered on (RA, Dec: 

J2000) = (12h56m26s,-05°49'30"), where no X-ray sources 
are found. Note that the background contamination is less 
than 1 % of source flux even in the faint X-ray states of the 
source. The data were rebinned to have at least 25 counts per 
bin, and the spectral fitting was performed using the energy 
range above 0.5 keV using XSPEC v.12.8.1. The Galactic col¬ 
umn density was fixed at 2.2 x 10^° cm“^. The data were an¬ 
alyzed and the flux and photon index were derived separately 
for each ObsID. A relation between the unabsorbed model 
flux (0.5-5 keV) and photon index is represented in Figure 5. 

Only the ObsIDs with an exposure of more than 600 s and 
with more than 7 spectral points (= d.o.f. > 5) were selected 
for the plot. A trend of a harder spectrum when the source is 
brighter is clearly detected. 

Data with ObsID of 80090001 and 35019132 were taken 
(quasi-) simultaneously with the NuSTAR observations, and 
here we report the details of those Swift-XRT observations. 

The observation with ObsID:80090001 was performed from 
2013 December 17 21:06:51 to 22:39:56 (UTC) with an ex¬ 
posure time of 2125 sec. Therefore, this observation did 
not exactly overlap with the NuSTAR observation, but it was 
the closest available, starting about 17 hours after the end 
of the first NuSTAR observation of the source. The spec¬ 
tral fit with a power-law model yielded a photon index of 
1.67 ±0.08 with a flux in the 0.5-5 keV band of (12.2 ± 

0.6) X 10“'^ erg cm“^ s“'. 

The other observation, with ObsID:35019132, was per¬ 
formed between 2014 January 01 00:20:26 and 22:50:54 
(UTC), which overlaps well with the second NuSTAR ob¬ 
servation. The exposure time of this Swift observation was 
6131s and the best-fit power-law model displayed a photon 
index of 1.42 ±0.03 with a flux in the 0.5-5 keV band of 
(19.1 ±0.3) X 10“'^erg cm“^ s“'. The spectral fitting results 
are reported in Table 3. Each snapshot observation during this 
ObsID was also analyzed separately. There were 15 snapshots 

http://heasarc.gsfc.nasa.gov/cgi-bin/W3Browse/swift.pi 
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Fig. 4.— X-ray light curves based on the count rates as measured by NuSTAR (black) and by Swiji-XRT (green). NuSTAR data are plotted in 1.5 hr bins, and 
the Swifi-XRT data are plotted for each snapshot. The vertical bars represent 1 cr statistical eiTors. 



Fig. 5.— Scatter plot of the X-ray flux (0.5-5 keV) vs. the X-ray photon 
index Fx of 3C 279 based on the Swifi-XRI data. The horizontal and vertical 
bars describe 1 a statistical errors for each axis. 

in total and the exposure time in each snapshot was about 
400 s typically, ranging from about 312 s to 594 s. The re¬ 
sultant count rates for all channels are shown in Figure 4. 

2.2.3. Joint spectral fit of NuSTAR and Swift-XRT data 

Joint spectral fits of the NuSTAR and Swift-XRT data were 
performed for each NuSTAR observation. As described in pre¬ 
vious sections, the data used for the spectral fitting were above 
0.5 keV for the Swift-XRT and 3-70 keV for NuSTAR. Here, 
we introduced a normalization factor of order (1 -3)% with 
respect to NuSTAR module A to account for differences in the 
absolute flux calibrations, and also to account for the offsets 
of the NuSTAR and Swift observing times (necessary for an 
analysis of a variable source). The Galactic column density 
was fixed at 2.2 x 10^° cm“^ as above. Simple power-law, 
a broken power-law and a double-broken power-law (for the 
second observation) models were used for the source spectral 
models. The joint spectral fitting results are summarized in 
Table 3. 

The joint spectrum during the first NuSTAR observation 
(Dec. 16) can be represented by a simple power law from 0.5 
to 70 keV with a photon index of 1.74 ±0.01 (x^/d.o.f. = 
688.6/679). The normalization factor of S'wi/f-XRT with re¬ 
spect to M/S'TAR module A is 1.01±0.05. The broken-power- 
law model improved the fit only marginally (~ 1 cr) with re- 



Fig. 6 .— Spectral energy distributions of 3C 279 in the soft-hard X-ray 
band based on the combined data from Swift-XRT and NuSTAR. The blue 
points show the results for observations on 2013 December 16-17 (in Period 
A), and the red points show the results for observations on 2013 December 31 
- 2014 January 01 (in Period C). The horizontal bars in data points describe 
the energy ranges of bins while the vertical bars represent 1 a statistical er¬ 
rors. 

sped to the power-law model. The result indicates that the 
X-ray spectrum of 3C 279 can be described by a single power 
law from the soft (0.5 keV) to the hard X-ray (70keV) band, 
which is supported by the results from the individual fits for 
each Swift-XRT and NuSTAR observation alone. 

For the joint spectrum during the second NuSTAR obser¬ 
vation, the simple power-law model did not result in an ac¬ 
ceptable fit, with x^/d.o.f. = 1072/996. Moreover, the nor¬ 
malization factor of Swift-XRT against the NuSTAR module 
A was ^ 0.7, which is clearly unacceptable. The broken 
power-law model, on the other hand, gave acceptable re¬ 
sults with x^/d.o.f. = 940.6/994, and the normalization of 
the 5w//r-XRT data was 0.97 ± 0.03. The break energy corre¬ 
sponded to 3.7 ±0.2keV with photon indices of 1.37 ±0.03 
and 1.76 ±0.01, respectively, below and above the break en¬ 
ergy. 

This break energy is located where the bandpasses of Swift- 
XRT and NuSTAR data overlap, corresponding respectively 
to the higher and the lower energy end of those data sets. 
We are confident that the spectral break is a real feature for 
the following reasons. There is a significant difference in the 
photon index in each Swift-XRT and NuSTAR data set con¬ 
sidered individually. This supports the conclusion that there 
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TABLE 3 

Parameters oe the spectral eits in X-ray band. 


Instrument 

(1) 

Tx/rxi 

(2) 

^brkl 

[keV] (3) 

rx2 

(4) 

^brkt 

[keV] (5) 

Tx3 

(6) 

const. 

XRT/module B (7) 

±2-I0keV 

(8) 

(9) 

F-test 

(prob.) 

Data on 2013 December 16 - 

17 (in Period A) 







Swift-XKl only 

1.67 ±0.08 





... / ... 

11.9 

21.10/18 (27.4%) 


NuSTAR only 

1.74 ±0.01 





■■■ /1.06±0.01 

11.0 

666.8/660 (41.9%) 


XRT 4- NuSTAR 

1.74±0.01 





1.01±0.05/1.06±0.01 

11.0 

688.6/679 (39.1%) 

- 


1 65+9.06 

4.5 ±0.7 

1.75 ±0.02 



l.llj;j[!^/1.06±0.01 

12.0 

686.2/677 (39.5%) 

30% (~ 1.0(7)“ 

Data on 2013 December 31 — 

2014 January 1 (in Period C) 







Swift-XRT only 

1.42 ±0.03 





... / ... 

23.7 

113.1/109 (37.5%) 


NuSTAR only 

1.75 ±0.01 





■■■ /1.00±0.01 

23.4 

831.2/886 (90.6%) 

- 


1.71 ±0.02 

9.3+]« 

1.81 ±0.02 



■■■ /1.00±0.01 

23.2 

820.8/884 (93.6%) 

0.4% (~ 2.9(7)“ 

XRT 4- NuSTAR 

1.73 ±0.01 





0.70±0.02/1.00±0.01 


1072.4/996 (4.6%) 



1.37 ±0.03 

3.7 ±0.2 

1.76±0.01 



0.97±0.03/1.00±0.01 

22.6 

940.6/994 (88.6%) 

- 


1 •1'7-KJ.U5 
*•^'-0.04 

-0.4 

1.72 ±0.02 

9.Mfl 

1.81 ±0.02 

0.97 ±0.03/1.00 ±0.01 

22.6 

933.4/992 (90.8%) 

0.22% (~ 3.1(7)'’ 

(log parabola)'^ 

1.39 ±0.03'' 

1 (fixed)® 

0.19±0.02f 



0.90 ±0.03/1.00 ±0.01 

22.7 

960.6/995 (77.8%) 



Note. — Col. (1): Instrument providing the data. Col. (2): Photon index for the power law model, or low-energy photon index for the broken power law model. Col. (3): Break energy 
[keV] for the broken power law model. Col. (4): High-energy photon index for the broken power law model. Col. (5): Second break energy. Col (6): Third index in the double-broken 
power law model. Col. (7): Constant factor of Swift-XKY/NuSTAR module-B data with respect to the NiiSTAR module-A data. Col. (8): Unabsorbed model flux in the 2-10 keV band, in 
units of 10“^^ [erg cm“^ s“*]. Col. (9): x^/degrees-of-freedom and a corresponding probability. 

a: compared to the simple-power-law model, b: compared to the broken-power-law model, c: logpar model in XSPEC. d: slope at the pivot energy, e: fixed pivot energy, f: curvature 
term. 


is a spectral break at an energy close to the overlap of the 
Swift-XRT and NuSTAR bandpasses. Each resultant individ¬ 
ual photon index is similar to the photon index derived from 
the joint fit below and above the break energy, respectively. 
The exposures of Swift-XRT and NuSTAR significantly over¬ 
lapped (see Figure 4), yielding a reasonable inter-calibration 
constant (0.97 ± 0.03). 

We also investigated a double-broken power-law model. 
The model yielded a probability of 0.22% (~ 3.1o-) that the 
improvement in the fit was due to chance against the bro¬ 
ken power-law model as assessed with an F-test. The second 
break energy appeared at 9.4^f4keV and the photon index 
became even softer above the second break energy, changing 
from 1.72 ±0.02 to 1.81 ±0.02. The spectral break at that 
energy was also seen in the fitting result for the NuSTAR data 
only at 9.3j;['3keV. All these results suggest that the X-ray 
spectrum during the high state (the second NuSTAR observa¬ 
tion) gradually softens with increasing energy, with the pho¬ 
ton index changing by ^ 0.4 from 0.5 keV to 70 keV. This 
is the first time that detailed, broad-band spectral X-ray mea¬ 
surements of 3C 279 show spectral softening with increas¬ 
ing energy. Furthermore the absence of spectral softening in 
the first (Dec. 16) NuSTAR observation clearly rules out the 
spectral shape as being caused by additional absorption. The 
joint X-ray spectral data points from the soft to the hard X- 
ray bands obtained by Swift-XRT and NuSTAR are plotted in 
Figure 6 in the E x F{E) (ergcm“^s“') form. Finally, a log 
parabola model was also tested using the logpar model in 
XSPEC. The pivot energy was fixed at 1 keV and the best-fit 
parameters are summarized in Table 3. The model also gave 
us acceptable fitting results, with x^/d.o.f. = 960.6/995. 

2.3. UV-Optical observations 
2.3.1. Swift-UVOT; UVbands 

The S’wi/f-Ultra-Violet/Optical Telescope (UVOT) data 
used in this paper included all observations performed dur¬ 
ing the time interval from 2013 November to 2014 April. The 
UVOT telescope cycled through each of the six optical and 
ultraviolet filters (W2, M2, Wl, U, B, V). The UVOT photo¬ 
metric system is described in Poole et al. (2008). Photometry 
was computed from a 5" source region around 3C 279 using 


the publicly available UVOT ftools data-reduction suite. 
The background region was taken from an annulus with inner 
and outer radii of 27."5 and 35", respectively. Galactic extinc¬ 
tion for each band in the direction of 3C 279 was adopted as 
given in Table 4. 

2.3.2. SMARTS: Optical-Near Infrared bands 

The source has been monitored for several years in the op¬ 
tical and NIR bands (B, V, R, J, and K bands) under the 
SMARTS project^^, organized by Yale University. Data re¬ 
duction and analysis are described in Bonninget al. (2012); 
Chatterjee et al. (2012), and the typical uncertainties for a 
bright source like 3C 279 are 1-2%. The publicly avail¬ 
able data were provided in magnitude scale. In a similar 
manner as presented in Nalewajko et al. (2012), the data in 
magnitude scale m\ were converted into flux densities as 
i^/v(ergs“'cm“^) = where Z\ is an effective 

zero point and A a is the extinction for each band. The ef¬ 
fective zero point is calculated as Z\ = 2.51og(c/Ae// x ff), 
where X^ff and f^ are parameters taken from Table A2 in 

TABLE 4 

Galactic extinctions in the UV-optical-Near IR bands as used in this 
PAPER. The extinctions are based on the reddening of 
E(B-V) = 0.029 MAG (SCHLEGEL ET AL. 1998) WITH Av/E(B-V) = 3.1. See 
ALSO Larionov et al. (2008). 


Band 

Ax 

Instruments 

W2 

0.271 

UVOT 

M2 

0.285 

UVOT 

Wl 

0.195 

UVOT 

U 

0.147 

UVOT 

B 

0.123 

UVOT, SMARTS 

V 

0.093 

UVOT, SMARTS, Kanata 

R,Rc 

0.075 

SMARTS, Kanata 

Ic, 

0.056 

Kanata 

J 

0.027 

SMARTS 

K 

0.010 

SMARTS 


http ; / /www . astro . yale . edu/smarts /glast / 
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Bessell et al. (1998). The extinctions were corrected using the 
values in Table 4. 

2.3.3. The Kanata telescope: Opticalphotopolarimetry 

We performed the V, Rq and /c-band photometry and Rq- 
band polarimetry observations of 3C 279 using HOWPol in¬ 
strument installed on the 1.5 m Kanata telescope located at 
the Higashi-Hiroshima Observatory, Japan (Kawabata et al. 
2008). We obtained 36 daily photometric measurements in 
each band, and 35 polarimetry measurements in the Rc band. 

A sequence of photopolarimetric observations consisted of 
successive exposures at four position angles of a half-wave 
plate; 0°, 45°, 22.5°, and 67.5°. The data were reduced us¬ 
ing standard procedures for CCD photometry. We performed 
aperture photometry using the APPHOT package in P YRAF 
and the differential photometry with a comparison star taken 
in the same frame of 3C 279. The comparison star is located 
at R.A. = 12:56:14.4 and Deck = -05:46:47.6 (J2000), and its 
magnitudes are T = 15.92, Rc= 15.35 (Bonning et al. 2012) 
and Ic= 14.743 (Zacharias et al. 2009). The data have been 
corrected for Galactic extinction as summarized in Table 4. 

Polarimetry with the HOWPol suffers from large instru¬ 
mental polarization ((5PD ^ 4%) caused by the reflection of 
the incident light on the tertiary mirror of the telescope. The 
instrumental polarization was modeled as a function of the 
declination of the object and the hour angle at the observation, 
and was subtracted from the observed value. We confirmed 
that the accuracy of instrumental polarization subtraction was 
better than 0.5% in the Rc band using unpolarized standard 
stars. The polarization angle is defined as usual (measured 
from north to east), based on calibrations with polarized stars, 
HD183143 and HD204827 (Schulz & Lenzen 1983). We also 
confirmed that the systematic error caused by instrumental po¬ 
larization was smaller than 2 ° using the polarized stars. 

2.4. Radio observations 
2.4.1. SMA: millimeter-wave band 

The 230 GHz flux density data was obtained at the Sub¬ 
millimeter Array (SMA), an 8 -element interferometer located 
near the summit of Mauna Kea (Hawaii). 3C 279 is included 
in an ongoing monitoring program at the SMA to determine 
the flux densities of compact extragalactic radio sources that 
can be used as calibrators at millimeter and sub-millimeter 
wavelengths (Gurwell et al. 2007). Observations of available 
potential calibrators are from time to time observed for 3 to 
5 minutes, and the measured source signal strength calibrated 
against known standards, typically solar system objects (Ti¬ 
tan, Uranus, Neptune, or Callisto). Data from this program 
are updated regularly and are available at the SMA website^^. 

3. MULTI-BAND OBSERVATIONAL RESULTS 
3.1. Light curve 

The multi-band light curves from the 7 -ray to the radio 
bands taken between MJD 56615 and 56775, are shown in 
Eigure 7 (covering the same period as in Eigure 1). The 7 -ray 
light curve measured by Ternu'-LAT is plotted using one-day 
time bins. The X-ray fluxes were measured by Swift-XRT in 

PYRAF is a product of the Space Telescope Sci¬ 
ence Institute, which is operated by AURA for NASA. See 
http : //www . stsci . edu/institute/software_hardware/pyraf 
http : / /smal . sma . hawaii . edu/cal list/cal list . html 


the 0.5-5 keV band. The third panel shows fluxes in the opti¬ 
cal V-band measured by 5wt/f-UVOT, SMARTS and Kanata 
as well as the R-band data measured by SMARTS and Kanata. 
The optical polarization data were measured by Kanata in the 
Rc-band. The 230 GHz fluxes were based on the results from 
SMA and also included results by ALMA^*. In the plot, the 
periods (A-D) as defined in Table 2 are also indicated. 

Generally, the source showed the most active states in the 7 - 
ray band at the beginning (including Period B, Elare 1) and the 
end (including Period D, Elare 3) of the epoch considered in 
this paper. In the X-ray band, we also see two high-flux states, 
in the first half and in the second half of this epoch. While 
in the first active phase the flux variation was not apparently 
well correlated between the 7 -ray and the X-ray bands, we 
can see flaring activities in both the 7 -ray the X-ray bands 
around Period D (^MJD 56750). 

During the epoch considered here, the optical flux showed 
significantly different behavior than that in the 7 -ray and the 
X-ray bands. In the beginning of this epoch, the measured 
fluxes were relatively low with relatively high polarization 
degrees, of ^ 20%. Around period B, the 7 -ray showed a 
very rapid flare with a hard photon index, but the source did 
not show any enhanced optical fluxes. After that, the optical 
fluxes started increasing gradually, with a drop of the polar¬ 
ization degree to ^ 5% after Period C. The 7 -ray and X-ray 
band fluxes dropped, but the optical flux still continued in¬ 
creasing, and peaked at simMJD 56720. In the largest flaring 
event in Period D, where the 7 -ray (> lOOMeV) and X-ray 
fluxes were highest, the optical flux showed only minor en¬ 
hancement, and had already started decreasing from its peak 
value. The optical polarization angle did not show any ro¬ 
tation throughout the observations considered here, and re¬ 
mained rather constant around 50° with respect to the jet di¬ 
rection observed by Very Long Baseline Interferometry ob¬ 
servations at radio bands (e.g., Jorstad et al. 2005). 

The 230 GHz flux was less variable compared to other 
bands, varying by about 50%, from ^ 8 Jy to 12 Jy. Even 
though the amplitude of the variation was much smaller, the 
general variability pattern of the 230 GHz band followed a 
similar pattern to that seen in the optical; a low state in the 
beginning of the epoch, followed by increased activity in the 
middle, and a decrease toward to the end of the interval. No 
prominent millimeter-wave flares corresponding to the large 
7 -ray flaring events (Elares 1-3) were observed. 

3.2. Spectral energy distributions 

Eigure 8 shows broadband SEDs for each period as de¬ 
fined in Table 2 (see also Eigure 7 in the light curves). The 
data sets include Remu'-LAT (see also Eigure 3), NuSTAR 
(for Periods A and C) and 5w//f-XRT (for Periods A, C, 
and D), Swift-VWOT (Wl for Period A, M2 and V for Pe¬ 
riod C, all six bands for Period D), SMARTS (B, V, R, J, 
K bands for all four periods) and Kanata {V, Rc bands for 
Period B), and SMA (for Periods A and C). All data in the 
figure were taken within the time spans as defined in Ta¬ 
ble 2 , but with a very slight offset in some optical data as fol¬ 
lows: for Period B (MJD 56646.4-56646.6), SMARTS data 
were taken during MJD 56646.348-56646.353 and Kanata 
data were taken starting at MJD 56646.8078. The SMARTS 
data observed during MJD 56750.1986-56750.2045 were for 
Period D (MJD 56750.210-56750.477). Unfortunately, no X- 

Taken from the ALMA Calibrator Catalogue, 
https : //almascience.eso.org/sc/ 
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ray observation was performed during Period B. For compar¬ 
ison, the SEDs during a polarization change associated with 
a 7 -ray flare observed in 2009 February, a low state in 2008 
August (Flayashida et al. 2012), and very-high-energy 7 -ray 
spectral points measured by MAGIC in 2006 (Albert et al. 
2008) are also included. 

4. DISCUSSION 

This campaign on blazar 3C 279 has three observational 
results of primary interest: ( 1 ) a rapid 7 -ray flare with very 
hard 7 -ray spectrum and no optical counterpart observed by 
Ferm/-LAT peaking at MID 56646 (Flare 1); (2) intraday vari¬ 
ability and a stable hard X-ray spectrum observed by NuSTAR 
in combination with significant X-ray spectral variations ob¬ 
served by Swift-XRT, and (3) a long trend of increasing op¬ 
tical flux without a corresponding increase in the 7 -ray flux. 
These results are discussed in detail in the following subsec¬ 
tions. 

4.1. Extreme ^-ray flare 

The 7 -ray flares peaking at MID 56646 (Flare 1, Period B) 
and MJD 56750 (Flare 3, Period D) are the brightest flares 
detected in 3C 279 by the Fermi-LAT. With photon fluxes 
of ~ 1.2 X 10“^ photons cm“^ s“* above 100 MeV, they are 
brighter by a factor ~ 4 than the flare peaking at MJD 54880 
(Abdo et al. 2010a), and comparable to the record fluxes de¬ 
tected by EGRET (Wehrle et al. 1998). 

Flare 1 was characterized by unprecedented rapid variabil¬ 
ity, with a flux-doubling time scale estimated conservatively 
at fvar — 2 h, and a very hard 7 -ray spectrum, with photon in¬ 
dex F-y ~ 1.7. Notably, we do not detect any simultaneous 
activity in the optical band. These facts make Flare 1 differ¬ 
ent from any previous 7 -ray activity observed in 3C 279. The 
closest analog was a rapid 7 -ray flare in PKS 1510-089 peak¬ 
ing at MJD 55854 with a flux-doubling time scale of ~ 1 h 
and photon index of F.^, ~ 2 (Saito et al. 2013; Nalewajko 
2013). However, in that case there were no simultaneous 
multiwavelength observations because of the proximity of 
PKS 1510-089 to the Sun. 

The SED for Flare 1 is presented in Figure 8 and 9 as Pe¬ 
riod B. The very hard 7 -ray spectrum measured by Fermi- 
LAT requires that the high-energy SED peak lies at energy 
> 2 GeV. We assume conservatively the actual SED peak lies 
at ~ 5 GeV, so that the apparent 7 -ray peak vF^ luminosity 
is L-y ~ 6 X 10"^* ergs“*. We further assume that this 7 -ray 
emission is produced by the external radiation Comptoniza- 
tion (ERG) mechanism (see Sikora et al. 2009 for a review 
of alternative mechanisms). The corresponding synchrotron 
component is expected to peak close to the optical band. The 
observed simultaneous optical/UV spectrum, and the lack of 
simultaneous optical variability, place very strong constraints 
on the Compton dominance parameter q = F^/L^yn > 300. 

In order to constrain the location along the jet r and the 
Lorentz factor Fj of the emitting region that produced Flare 
1, we use the recent model of Nalewajko et al. (2014). In that 
model, the allowed parameter space for the 7 -ray emitting re¬ 
gion is defined by three constraints: ( 1 ) a jet collimation con¬ 
straint FjSj < 1 , where is the jet half-opening angle, ( 2 ) a 
constraint on the SSC luminosity Lssc < Tx, where Lx is the 
observed X-ray, hard X-ray or soft 7 -ray luminosity (depend¬ 
ing on the expected energy of the SSC peak), and (3) a con¬ 
straint on the radiative cooling time scale Lcooi < lOOMeV, 
where Lcooi is the characteristic observed energy of the 7 - 
ray photons produced by the electrons for which the radiative 


cooling time scale is comparable to the variability time scale. 
The proper size of the emission region is estimated directly 
from the observed variability time scale R ~ VcRar .obs/(l 4“ fl), 
and it is related to the jet opening angle by R/r = 9y In ad¬ 
dition to the parameters discussed above, we adopt a stan¬ 
dard ratio of the Doppler-to-Lorentz factors D/Fj = 1, and 
we also need to specify the upper limit on the expected syn¬ 
chrotron self-Compton (SSC) component. As the SSC com¬ 
ponent likely peaks in the MeV band, where we do not have 
any observational constraints, we will conservatively assume 
that Lssc < L-y/30 ~ 2 x 10"^^ ergs“^ The constraints on the 
parameter space are shown in Figure 10. The yellow-shaded 
area indicates the allowed location of the 7 -ray emitting re¬ 
gion. For reasonable values of the jet Lorentz factor 20 < 
Fj < 30, it should be located between 0.015 pc < r < 0.15 pc, 
where the external radiation is dominated by the broad emis¬ 
sion lines (cblr — 0.8 pc). At larger distances, the energy den¬ 
sity of external radiation fields will be insufficient to provide 
efficient cooling of the electrons producing the 100 MeV pho¬ 
tons on the observed time scales, and also it will be more dif¬ 
ficult to maintain a sufficiently high energy density to power 
such a luminous 7 -ray flare from a very compact emission 
region. The Lorentz factor should be Fj > 17, although this 
constraint would be stronger if we assumed a lower Lssc- For 
the subsequent modeling of the SED, we will adopt two pos¬ 
sible solutions indicated in Figure 10: (Bl) r = 0.03 pc and 
Fj = 20, and (B2) r= 0.12 pc and Fj = 30. This corresponds to 
the jet opening angle 0j satisfying: (Bl) Fj^j = 0.61 and (B2) 
rj6lj=0.34. 

We use the Blazar code (Moderski et al. 2003) to model 
the SED of 3C 279 during Elare 1 with a standard leptonic 
model including the synchrotron, SSC and ERC processes. 
The distribution of external radiation is scaled to the accretion 
disk luminosity of Ld ~ 6 x 10 "^^ ergs“' using standard rela¬ 
tions for the characteristic radii of the broad-line region (BLR) 
and the dusty torus (Sikora et al. 2009) with covering factors 
■Cblr = Cir = O-F In general, we use a double-broken power- 
law energy distribution of injected electrons Niy) oc 7 “^' with 
two breaks at 71 and 72 and three indices: p\ for 7 < 71 , pi for 
7 i < 7 < 72 , and pj, for 7 > 72 . In order to reproduce a very 
hard 7 -ray spectrum in the fast-cooling regime, we set pi = 1. 
The goal of the modeling is to match the ERC(BLR) peak 
with the observed 7 -ray peak by adjusting the maximum elec¬ 
tron Lorentz factor, 71 , and the electron distribution normal¬ 
ization, and then to match the synchrotron component with the 
optical/UV data by adjusting the co-moving magnetic field B'. 
One should remember that because of the lack of any optical 
activity simultaneous with Elare 1, the synchrotron compo¬ 
nent should actually be below the optical/UV data points, so 
B' should be treated more like an upper limit than an actual 
value. The results of the SED modeling are presented in Fig¬ 
ure 9, and essential model parameters are listed in Table 5. 
For solution Bl we obtain 71 = 3700 and B' = 0.31 G, and for 
solution B2 we obtain 71 = 2800 and B' = 0.3 G. 

We consider the basic energetic requirements for produc¬ 
ing such an SED. We can estimate the total required jet power 
as Lj ~ L..y/( 77 jr?), where pj ~ 0.1 is the radiative efficiency 
of the jet. And we can use the estimated magnetic field 
strength to calculate the magnetic jet power Lb = ttR^V'^u'^c = 
(c/8)(rj0j)^(BV)^, where R = 9/ is the jet radius, is the 
jet half-opening angle, and i/g = B'^/( 87 r) is the magnetic en¬ 
ergy density. For solution B1 we obtain Lj ~ 1.5 X 10"^^ ergs ^ 
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Fig. 7.— Multiwavelength light curves of 3C 279 covering the same period as in Figure 1. From the top, the panels show: (1) 7 -ray photon flux above 100MeV 
in 1-day bins from Fermi-LP[Y\ (2) X-ray flux density between 0.5-5 keV from Swift-XKI; (3) optical flux density from 5w(yr-UVOT, SMARTS (Ki?), and Kanata 
(k"); (4) optical polarization degree (scale on the left) and electric vector polarization angle (EVPA, scale on the right) from Kanata; (5) mm flux density (230 GHz) 
measured by SMA and ALMA. The vertical dotted lines indicate the periods (A—D) as defined in Table 2 when the 7 -ray spectra were extracted. A gap in the 
7 -ray data by Fermi-hA^ around MJD 56680-56690 is due to a ToO observation of the Crab Nebula, during which time no exposure was available in the direction 
of 3C 279. The vertical bars in data points represent 1 a statistical errors and the down arrows indicate 95% confidence level upper limits. 


and Lb — 1-1 X 10"^^ ergs“', and for solution B2 we obtain 
Lj ~ 7 X 10"^^ ergs“' and Lb — 5 x 10"^^ ergs“'. In both cases, 
the required magnetic jet power is a tiny fraction of the to¬ 
tal jet power. This fraction is higher for solution B2, with 
Lb/Lj ~ 0.7 X 10“^. This indicates that the emitting region 
responsible for Flare 1 is very strongly matter-dominated 
(cf. Janiaketal. 2015), although several observed SEDs of 
3C 279 in 2008-2010 can be described by an equipartition 
model (Dermer et al. 2014). 

The total required jet power Lj appears very high compared 
to the accretion disk luminosity La. For solution B1 we ob¬ 
tain Lj/Ld ~ 25, and for solution B2 we obtain Lj/La ~ 11. 
We note that there is no signature of increased disk luminos¬ 
ity in the UV data even for the lowest-flux state (Period A), 


which gives La < 9 x 10“*^ ergs“^ Assuming the black hole 
mass of Mbh ~ 5 x IO^Mq (see Section 1), the Eddington lu¬ 
minosity is LEdd — 8 X 10^® ergs“*, hence La/Lead — 0.08 and 
Lj/Ledd — 0.9 for solution B2. Solution B1 with moderate jet 
Lorentz factor Ej = 20 predicts a super-Eddington jet power. 
The predicted jet power can also be decreased when assum¬ 
ing a higher jet radiative efficiency ryj >0.1. Taking this into 
account, the energetic requirements for this flare are consis¬ 
tent with the typical relation between Lj and La for FSRQs 
(Ghisellini et al. 2014). 

In the standard scenario of magnetic acceleration of rel¬ 
ativistic black-hole jets, the jets are initially dominated by 
the magnetic energy, which is gradually converted into the 
kinetic energy of plasma until the jet magnetization a ~ 
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Fig. 8 .— Broadband spectral energy distributions of 3C279 for the four observational periods defined in Section 2.1 (see also Table 2 and Figure 7). The 
vertical bars in data points represent 1 a statistical eiTors and the down an'ows indicate 95% confidence level upper limits. The plot includes historical SEDs of 
3C279 in a low state (in 2008 August) and in a flaring state (in 2009 February) from the 2008-2010 campaign (Hayashida et al. 2012). The measured spectral 
fluxes by MAGIC in 2006 are also plotted (Albert et al. 2008). 


Ls/iLj-LB) ^ 1 (Begelman & Li 1994; Komissarov et al. 
2007). Additional conversion of the jet magnetic en¬ 
ergy is possible locally under special conditions, e.g. rar¬ 
efaction acceleration induced by rapid decline in external 
pressure (Tchekhovskoy et al. 2010; Komissarov et al. 2010; 
Sapountzis & Vlahakis 2014). However, this additional jet 
acceleration is likely to lead to the loss of causal contact 
across the jet, and very wide opening angles (Lj 0j> 1), more 
typical of gamma-ray bursts. Magnetic energy can be dissi¬ 
pated directly in the process of magnetic reconnection. How¬ 
ever, it is unlikely that relativistic reconnection can operate 
with sufficient efficiency to convert highly magnetized plasma 
with tr ^ 1 to very weakly magnetized plasma with a ^ 10“"^. 
Additional matter could be injected into the jet by individ¬ 
ual massive stars crossing the jet (Khangulyan et al. 2013), 
or at the jet base during a brief pause in the jet production 
(Dexter et al. 2014). The latter possibility potentially allows 
for a more uniform distribution of an unmagnetized plasma 
layer across the jet. 

The very hard electron energy distribution with pi ~ 1 ex¬ 
tending up to 7 > 2000 , required to explain the very hard 
7 -ray spectrum of Flare 1, is challenging for many par¬ 
ticle acceleration mechanism and emission scenarios (e.g., 
Blandford & Levinson 1995). Very hard electron spectra 
can be obtained in relativistic magnetic reconnection, but 
they require extremely high electron magnetization CTc > 100 
(Sironi & Spitkovsky 2014; Guo et al. 2014; Werner etal. 
2014). In the case of an electron-proton jet composition with 


no e'^e~ pairs, one has cr^ ~ (^p/^e){inp/ine)a, where jp and 
7 e are typical Lorentz factors of protons and electrons, respec¬ 
tively. In principle, it is possible that <7^ ^ cr, so that such ex¬ 
treme acceleration of electrons is possible even in the case of 
(7 < 1. The final outcome of the acceleration depends on how 
the dissipated magnetic energy is shared between the protons 
and electrons; the first study of relativistic electron-ion recon¬ 
nection suggests roughly equal energy division (Melzani et al. 
2014). 

We propose that Flare 1 of 3C 279, together with the simi¬ 
lar flare of PKS 1510-089 peaking at MID 55854 (Saito et al. 
2013), constitute an emerging class of rapid 7 -ray events char¬ 
acterized by flux-doubling time scales of a couple of hours, 
very hard 7 -ray spectra with spectral peaks in the GeV band, 
and significant time asymmetry with longer decay time scales 
(Nalewajko 2013). Moreover, the results of this work indi¬ 
cate that such events do not have significant multiwavelength 
counterparts. Since only two clear examples were identified 
in bright blazars during ^ 6 years so far of the Fermi mission, 
they appear to be rare events, and may not represent typical 
conditions of dissipation and particle acceleration in blazar 
jets. 

In Figure 9, we also present two SED models for Flare 3 
(Period D). This flare is characterized by a typical 7 -ray spec¬ 
trum, and a more typical Compton dominance, as compared 
to Flare 1. In addition, for Flare 3 we have simultaneous 
UV and X-ray data from Swift. The soft X-ray spectrum is 
very hard, with Fx = 1.22 ±0.07. We first attempted — in 
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Fig. 9.— Left panel. Spectral energy distributions of 3C 279 during the brightest 7 -ray flares — Flare 1 (Period B, red points) and Flare 3 (Period D, blue 
points) — see Section 4.1 for discussion. Right panel. SEDs during two NiiSTAR pointings — Period A (orange points) and Period C (magenta points) — see 
Section 4.2 for discussion. Solid and dashed lines show SED models obtained with the leptonic code Blazar. Model parameters are listed in Table 5. Black and 
gray lines show historical data and SED models from Hayashida et al. (2012). Black dashed line shows the composite SED for radio-loud quasars (Elvis et al. 
1994) normalized to Ld = 6 x lO"^^ ergs“^. The inset illustrates schematically the decomposition of each SED model into contributions from individual radiative 
mechanisms: (in order of increasing peak frequency) synchrotron, SSC, ERC(IR), ERC(BLR). The axes and line types are the same as in the main plot. 
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Fig. 10.— Constraints on the parameter space of location r and Lorentz factor Fj for the emitting region producing 7 -ray Flare 1 (Period B) (see Nalewajko et al. 
2014, for detailed description of the model). The following constraints are shown: from jet collimation parameter Fjfij (solid red lines), from SSC luminosity 
Lssc (dashed blue lines), from energy threshold for efficient cooling £cool (dotted magenta line), from the characteristic wavelength of synchrotron self-absorption 
Assa (dot-dashed orange line), from the minimum required jet power Lj (double-dot-dashed green lines), and from the characteristic energy of intrinsic 7 -ray 
absorption iJmax.obs (solid gray line). The region allowed by the collimation constraint (Fj^j < 1), the SSC constraint (Lssc < 2 X lO'^^ergs^*) and the cooling 
constraint (Lcool < 100 MeV) is shaded in yellow. Two particular' solutions B1 and B2 for which SED models are shown in Figure 9 are indicated. 


model D1 located in the BLR — to explain this X-ray spec¬ 
trum by SSC emission from a very hard electron energy dis¬ 
tribution (pi = 1). By coincidence, model B1 described in 
the previous subsection does exactly that. However, since the 
7 -ray spectrum for Period D is much softer than the excep¬ 
tionally hard 7 -ray spectrum for Period B, in model D1 we 
need to adopt a break in the electron energy distribution at 
7 br ^ 200, which shifts the observed peak of the SSC compo¬ 
nent to ^ 100 keV; hence the X-ray part of the SSC spectrum 
is too soft to explain the observed X-ray spectrum. We note 
that Paliyaetal. (2015) present an SED model for a period 
overlapping with our Period D, in which the X-ray spectrum 
is matched with the SSC component. They made the model 


by adopting a higher value of 7 br, which requires a superpo¬ 
sition of ERC(BLR) and ERC(IR) components to explain the 
7 -ray spectrum. In model D2 we were able to explain the X- 
ray spectrum with the low-energy tail of the ERC emission. 
The emission region in model D2 is located outside the BLR, 
and the entire high-energy component is strongly dominated 
by the ERC(IR) emission. We adopted a higher jet Lorentz 
factor Lj = 30 and the low-energy electron distribution index 
Pi = 1.6 for 7 < 100 (see Table 5). While model D2 matches 
the observed X-ray spectrum, because it is located at relatively 
large distance r ~ 1.1 pc, it predicts a rather long observed 
variability time scale of fvar,obs 2 d. 

Such extremely high flux spectra as observed in Periods B 
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TABLE 5 

Parameters of the SED models presented in Eig. 9. 


Model 

A 

B1 

B2 

C 

D1 

D2 

r [pc] 

1.1 

0.03 

0.12 

1.1 

0.03 

1.1 

Tj 

8.5 

20 

30 

10.5 

25 

30 

Ej^j 

1 

0.61 

0.34 

1 

1 

1 

B' [G] 

0.13 

0.31 

0.3 

0.13 

1.75 

0.14 

Pi 

1 

1 

1 

1 

1 

1.6 

71 

1000 

3700 

2800 

1000 

200 

100 

pi 

2.4 

7 

7 

2.4 

2.5 

2.5 

72 

3000 

— 

— 

3000 

2000 

6000 

P3 

3.5 

— 

— 

3.5 

5 
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and D could in some cases extend to even higher energies, 
and may possibly be detectable by ground-based Cherenkov 
telescopes (MAGIC, H.E.S.S., VERITAS, CTA). Despite its 
moderate redshift, 3C 279 was detected twice by MAGIC 
(Albert et al. 2008; Aleksic et al. 2011a) before the Fermi era. 
In this work, we have argued that the 7 -rays originate at a ra¬ 
dius ^ 0.1 pc, which is comparable with the estimated size 
of the broad emission line region tblr based upon reverber¬ 
ation mapping campaigns of other AGN (e.g., Bentz et al. 
2006; Kaspi et al. 2007). This radius is also roughly compa¬ 
rable to the minimum radius from which the highest energy 
photon observed during our campaign — Eobs = 26.1 GeV 
{E ^ 40 GeV in the quasar rest frame) in Period A — can 
escape without pair production. The highest energy photons 
detected in Periods B and D were 10.4 GeV and 13.5 GeV, 
respectively (see Eigure 1). We could not distinguish if the 
non-detection of > 15 GeV photons was due to the absorption 
by the BLR photons or just the statistical limitation of the 
short integration time for the spectra. Our emission models 
indicate a sharp drop in the source intrinsic spectral shape at 
>10 GeV energies due to adopting a very steep high-energy 
electron distribution index p 2 = 7. In addition, the 7 -ray emis¬ 
sion above 10 GeV produced in the ERC(BLR) process is sup¬ 
pressed due to the reduction of the scattering cross section in 
the Klein-Nishina regime. 

The importance of 7 -ray absorption in the pair produc¬ 
tion process depends on the abundance of soft photons pro¬ 
duced in the jet environment. One source of soft photons 
is the emission lines radiated by the broad emission line 
clouds. The optical depth depends on the geometrical shape 
of the BLR, and is significantly reduced for the flat geome¬ 
tries (Tavecchio &. Ghisellini 2012; Stern &. Poutanen 2014). 
Specifically, the results of Tavecchio & Ghisellini (2012) cal¬ 
culated for Ld = 5 X 10^^ ergs“* and for an intermediate geo¬ 
metric case indicate that absorption from BLR photons is not 
significant for 7 -ray photons with Eobs 50 GeV emitted at 
r ~ rBLR, and those with Eobs ^ 20 GeV at r ^ O.Itblr- How¬ 
ever, other models of the BLR can be expected to set a larger 
lower bound on the emission distance scale r. Another source 
of soft photons is Thomson scattering by the hot inter-cloud 
medium that is commonly invoked to confine the clouds by 
ram pressure. If we adopt the spectral component associated 
with the accretion disk by Pian et al. (1999), then the radius 
of the ‘ 7 -sphere’ (for the 40 GeV photons in the source 
frame) is Q.2{{tt) /0.01)pc, where (tt) is the mean Thom¬ 
son depth (e.g., Blandford & Levinson 1995). 

The difficulties faced by the leptonic models do not neces¬ 
sarily mean that they should be abandoned, unless there ex¬ 


ists a better alternative. Hadronic models have been applied 
to 3C 279 (Bottcher et al. 2009; Petropoulou & Mastichiadis 
2012; Diltz et al. 2015), however, they always require an ex¬ 
tremely large jet power of order 10 "^® ergs“*, which is difficult 
to reconcile with observations of radio galaxies and theories 
of jet launching (Zdziarski & Boettcher 2015). 

4.2. Spectral variability in X-ray and hard X-ray bands 

The X-ray and hard X-ray behavior of 3C 279 revealed by 
Swift-XRT and NuSTAR appears to be quite complex. This 
is especially interesting since the origin of this emission in 
LSRQ-type blazars is in general controversial and poorly un¬ 
derstood (Sikoraet al. 2013). Here, we focus on constraining 
the mechanism of X-ray and hard X-ray emission observed 
during the two NuSTAR pointings (Periods A and C). 

The X-ray flux observed during the NuSTAR pointings is 
relatively high for 3C 279, higher by a factor of ~ 5 than typ¬ 
ical X-ray fluxes measured during the 2008-2010 campaign 
(Hayashida et al. 2012). In Period C, the X-ray flux is higher 
by a factor of ^ 2 than in Period A, and it also shows some 
intraday variations. A comparison of the broadband SEDs for 
Period A and C suggests a roughly linear relation between 
optical. X-ray and 7 -ray variations. However, this may be 
misleading, as the optical flux shows a long-term systematic 
increase that is not evident in the X-ray and 7 -ray light curves. 
With typical photon index Lx — 1.7, the observed X-ray spec¬ 
trum is relatively hard, with a spectral break observed in Pe¬ 
riod C at ~ 3.5 keV. The Swift-XRT data show a clear anti¬ 
correlation between the soft X-ray photon index and the soft 
X-ray flux (see Ligure 5). 

The X-ray emission of LSRQs is typically attributed either 
to the SSC emission of medium-energy electrons, or to the 
ERG emission of low-energy electrons. Adopting a one-zone 
model in the SSC scenario for the X-ray band, we attempted 
to explain the observed broadband SEDs for Periods A and 
C together with the synchrotron emission for the optical band 
and the ERC component for the 7 -ray band. Eigure 9 shows 
the SEDs for Periods A and C with the emission models based 
on the parameters in Table 5. Actually, the observed high X- 
ray flux and hard X-ray spectra challenge the SSC scenario. In 
order to match the relatively high observed X-ray luminosity 
with that expected to be produced via SSC, as well as the cor¬ 
responding synchrotron component and the ERC components 
with a condition of Ej = 1 , one needs to adopt a rather low jet 
Lorentz factor Ej ~ 8.5 for Period A and Ej ~ 10.5 for Period 
C. A low jet Lorentz factor requires a more powerful jet; in the 
case of Period C with 7 -ray luminosity oe. 4 x 10^^^ ergs”' 
we estimate Ej ~ E-y/(pjE?) ~ 4 x 10^®(77j/0.1)”' ergs”*. The 
jet Lorentz factor can be higher, and the required jet power 
lower, by allowing that Ej0j < 1. Typical X-ray SEDs from 
the SSC component in ESRQs are flat; it is possible to ob¬ 
tain a hard SSC spectrum by choosing an electron energy 
distribution peaking at 7 peak ^ 300. However, in order to re¬ 
produce apparent linear relations between optical. X-ray and 
7 -ray fluxes from Periods A to C, it was necessary to ad¬ 
just the value of Ej to compensate for the quadratic depen¬ 
dence of the SSC luminosity on the synchrotron luminosity 
Esse C!C Ejyjj. Therefore, while it is possible to make a one- 
zone model that fits the observed SEDs for both Periods A 
and C with X-rays produced by the SSC process and with rea¬ 
sonable jet power, such a model cannot naturally account for 
the observed flux variations over the multiwavelength bands. 
We note that lower Lorentz factors Ej ~ 10 required for mod- 
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eling the SEDs for Periods A and C, together higher Lorentz 
factors Ej ~ 30 required for modeling the SEDs for Periods 
B and D, could indicate the existence of a spine-sheath jet 
structure (Ghisellini et al. 2005). In such case the 7 -ray flares 
would be produced in the fast spine and the bulk of X-ray and 
hard X-ray emission in the slow sheath. 

In the ERC scenario, on the other hand, the observed X-ray 
flux could be dominated by the low-energy tail of either the 
ERC(BLR) or ERC(IR) components, depending on the loca¬ 
tion of the emitting region. In one-zone models, this would 
require the X-ray spectra to be related to the 7 -ray spectra by 
a single spectral component. This would explain the appar¬ 
ent linear relation between the X-ray and 7 -ray data, although 
there are insufficient X-ray observations in the current cam¬ 
paign to probe the direct correlation between X-ray and 7 -ray 
fluxes. (During the previous campaign on 3C 279, no signif¬ 
icant correlation was detected between the X-ray and 7 -ray 
fluxes; Hayashidaet al. 2012). Judging from the simultane¬ 
ous X-ray, hard X-ray and 7 -ray SEDs, it could be possible to 
connect them by a single spectral component, especially for 
Period A. This would require a power-law extension of the X- 
ray spectrum all the way to the low-energy end of the 7 -ray 
spectrum. This is possible only for the ERC(IR) component 
where no cooling break is expected. However, because the 
low-energy ERC(IR) emission would be produced deep in the 
slow-cooling regime, very little flux variability would be ex¬ 
pected on daily time scales. The ERC(BLR) component is 
very likely to feature a cooling break, and possibly an addi¬ 
tional low-energy spectral break at 2.6(rj/20)^ keV produced 
by trans-relativistic electrons (7 ^ 1 ). 

We conclude that there is no one-zone leptonic SED model 
that can satisfactorily explain the production of X-ray emis¬ 
sion observed by NuSTAR in 3C 279. Various alternative 
mechanisms can be proposed where the observed X-ray emis¬ 
sion is produced at a different location from the optical and/or 
7 -ray emission. Observation of a transient spectral break 
at ~ 3.5 keV in Period C may indicate a superposition of 
two spectral components in the X-ray band. Similar spec¬ 
tral breaks observed in high-redshift ESRQs were interpreted 
as due to very strong absorption (Eabianet al. 2001). How¬ 
ever in our case this interpretation is challenged by the lack 
of a break in Period A, only two weeks earlier. Alterna¬ 
tive mechanisms for the X-ray emission include IC emission 
from the accretion-disk corona, ERC or synchrotron emission 
from the jet acceleration region, and hadronic mechanisms 
(Bottcher et al. 2009). However, for most mechanisms it may 
be challenging to explain the relatively high X-ray luminos¬ 
ity Lx — 2 X 10"^® ergs“' observed during Period C, a factor 
~ 3 higher than L^. This certainly excludes the accretion-disk 
coronal emission proposed tentatively for PKS 1510-089 in 
the low state (Nalewajko et al. 2012). The ‘jet base’ scenario, 
where luminous weakly beamed X-ray emission is produced 
at very short distances from the supermassive black hole at 
which the jet is not yet fully accelerated, is motivated by re¬ 
cent observations of misaligned ESRQs (Bostrom et al. 2014), 
and will be investigated in detail elsewhere. 

4.3. Optical behavior 

The observed optical flux shows a systematic increase by a 
factor of ^ 4 over a period of ^ 80 days (MJD 56645-56725). 
Superposed on this trend are weak flares that correspond very 
poorly to the strong 7 -ray flares. No similar systematic trend 
is seen in the X-ray and 7 -ray light curves. This is in con¬ 


trast with the good overall correlation between optical and 7 - 
ray fluxes in 2008-2010 (Hayashida et al. 2012). The lack of 
overall correlation between the optical light curve and the ra¬ 
dio (mm-band). X-ray, and 7 -ray light curves suggests that 
they are produced by different populations of electrons, and 
most likely at different locations. Moreover, the apparent lin¬ 
ear relation between SEDs for Periods A and C is merely a 
coincidence. As the overall radiative output from 3C 279 is 
always dominated by the 7 -ray emission, the systematic long¬ 
term increase in the optical flux could be due to a system¬ 
atic increase of the magnetic held strength, or a systematic 
decrease of the external radiation energy density. The for¬ 
mer option would be problematic if the jet magnetization re¬ 
mains constant, which would lead to an increase in the total jet 
power, and ultimately to an increase in the 7 -ray luminosity, 
which is not observed. 

5. CONCLUSIONS 

We report the results of observations of the well-studied 
7 -ray luminous blazar 3C 279 at the end of 2013 and be¬ 
ginning of 2014, when the object entered a bright and active 
state. The Lermi-LAT observations revealed multiple distinct, 
bright flares, and recorded the highest 7 -ray flux state of the 
source since the launch of Fermi, at F{E > lOOMeV) reaching 
10“^ photons cm“^ s“^ on 2013 December 20 and 2014 April 
03. The high flux of the source allowed us to integrate the 
7 -ray flux on time scales as short as one Lermi-LAT orbit (96 
min). This in turn allowed us to establish the variability time 
scales to be as short as Lar — 2 hours. One of these flares re¬ 
vealed an unprecedentedly hard 7 -ray spectrum, with a pho¬ 
ton index E.^ ~ 1.7, unusual for this kind of source. 

Two NuSTAR observations provided the first precise mea¬ 
surement of the hard X-ray spectrum of 3C 279 up to 
70 keV. The NuSTAR observations were complemented by 
more frequent Swift observations. The best-fit model for 
the joint spectra by NuSTAR and Swift-XRT during the first 
NuSTAR observation (2013 December 16) was consistent 
with a simple power law, which usually has been observed 
in past X-ray observations involving Suzaku and XMM- 
Newton (Hayashida et al. 2012). On the other hand, the sec¬ 
ond simultaneous NuSTAR and Swift observations at the end 
of 2013 revealed a spectral structure that was harder (Exi — 
1.37) below ~ 3.5 keV and softer (Ex 2 — 1.76) above that 
energy. Such a spectral structure has not been observed in 
3C 279 in any previous X-ray observations. In addition, the 
second NuSTAR observation (2013 December 31) indicated 
an increase of the X-ray flux by ~ 50% during the 1-day 
pointing. In the soft X-ray data from Swift-XRT, we found 
a clear correlation between the flux and photon index with 
a harder-when-brighter trend. More detailed studies offer po¬ 
tential for better understanding of the origin of X-ray emission 
in ESRQ blazars. 

The optical flux of the source steadily increased since the 
beginning of this multiwavelength campaign, but did not show 
clear, large-amplitude flares such as those seen in 7 -rays. The 
optical flux of the source does not appear clearly correlated 
with the 7 -ray flux, in contrast to the behavior measured in 
2008-2010 (Hayashida et al. 2012). It is possible that the op¬ 
tical flux might be arising in multiple locations along the jet, 
as expected if the electrons radiating in the optical are of dif¬ 
ferent energy than those making the X- and 7 -rays, which is 
suggested by the apparent lack of correlation between the op¬ 
tical flux and the optical polarization degree. 

We modeled the broad-band spectrum of the source at mul- 
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tiple epochs, including two very bright 7 -ray flares, and the 
two NuSTAR pointings. The very hard 7 -ray spectrum and 
very high Compton dominance during the first 7 -ray flare 
are very challenging to explain in the standard one-zone syn¬ 
chrotron plus SSC/ERC model, requiring electron energy dis¬ 
tribution index p ~ 1 , a high jet power compared to the ac¬ 
cretion disk luminosity Lj/La ^ 10 , and a very low magnetic 
fraction of the jet power Lb /Lj < 10“^. In addition, no single¬ 
zone modeling of any single epoch of the broad-band SED 
can satisfactorily explain the production of X-ray emission at 
the same time as it explains the optical and 7 -ray emission. 
This conclusion is consistent with the finding that variations 
in the X-ray flux of 3C 279 are not always simultaneous with 
variations in the 7 -ray or optical fluxes. 
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